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ABSTRACT 

We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (bands 6b) and 1788.4 ~ 
1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines 
found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread 
emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines 
is significantly diminished in these bands compared to lower frequency HerschelfHIFl bands. 
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^ 1. Introduction 

■ The Kleinmann-Low nebula within the Orion molecular cloud 
(Orion KL) is the best studied massive star forming region in 

. the Mil ky Way. This re gion is characterized by a high IR lumi- 
' nosity (iKleinmann & L ow 1967) and rich molecular line emis- 
sion. As such, it has been the subject of numerous molecular line 
surveys in the millimeter and submillimeter that have character- 
' ized its mm/sub-mm wave spectrum (see e.g. lSchilke et al.ll997t 
, IComito et al.ll2005t iTercero et al.ll201 0', and references therein). 
' These surveys reveal the presence of a prodigious variety of 
molecular species in addition to several distinct spatial/velocity 
1 components ( i.e. the hot core, compact ridge , plateau, and ex- 

■ tended ridge; IB lake etal.|[l 9871: iPersson et al] l2007). These ob- 
' servations provide insight into the complex chemical and physi- 
cal processes that characterize how massive stars form and in- 
teract with their natal environment. Because Orio n KL is the 

, close st such massive star forming region (~414 pc; Menten et al.l 
l2007h . it is an ideal choice for further inquiry in understanding 
the chemistry and physics of the gas in close proximity to these 
stars. 

Although the subject of much spectroscopic study in the mm 
/sub-mm [A > 300;um) during the past 30 years, high resolution 
observations at Terahertz (THz) frequencies of Orion KL have 
been unavailable from the ground due to atmospheric absorp- 
tion. ISO provided the first comp rehensive spectrosc opic view 
of Orion KL at these wavelengths. iLerate et al.l(l2006h presented 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



spectroscopic observations in the wavelength range 44 - 188 
//m (1.6 - 6.8 THz) with a resolving power of A/AA ~ 6800- 
9700 using the long wavelength spectrometer (LWS) on board 
ISO. These data showed a spectrum dominated by emission from 
H2O, OH, and CO, but little or no emission from more complex 
species (such as methanol, methyl formate, dimethyl ether, etc.), 
which litter the spectrum at submillimeter wavele ngths. 

The HIFI instrument (Ide Graauw et Si\!'201(t) on board the 
Herschel Space Observatory (Pilbratt et al. 2010) provides the 
first opportunity to characterize the THz spectrum with high 
spectral resolution and sensitivity. In this Letter, we present 
the first high resolution (A/AA ~ 10^) spectrum of Orion KL 
above 1 .57 THz obtained using the HIFI instrument. These ob- 
servations, taken as part of the guaranteed time key program 
Herschel observations of EXtra-Ordinary Sources: The Orion 
and Sagittarius B2 Starforming Regions (HEXOS), are able to 
probe the chemical inventory and kinematic structure of Orion 
KL at an unprecedented level. In this work , we characterize the 
high resolution THz spectrum. We further demonstrate and dis- 
cuss why the observed Une density is reduced when compared to 
lower frequencies. 



2. Observations 

The observations were carried out on March 22-23, 2010 using 
the wide band spectrometer (WBS) with a spectral resolution of 
1.1 MHz (0.19 km s' at 1.75 THz) over a -2.4 GHz IF band- 
width. The data were taken in dual beam switch (DBS) mode 
using the fast chop setting pointed towards the Orion hot core at 
coordinates ffyjooo = 5''35"'14.5' and 6j2(m = -5°22'30.9". The 
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Fig. 1. SSB spectrum of the Orion KL hot core in Bands 6b (top panel) and 7b (bottom panel) smoothed to a velocity resolution of 
~4.5 km s"' . The strongest lines (Tmb ^ 7 K) are labelled. 



beam size at 1.75 THz is 12" and the DBS reference beams lie 
approximately 3' east and west. Both H and V polarization data 
were obtained. However, we only present the H polarization here 
because the mixer is optimi zed for th e H polarization. These data 
were reduced using HIPE (iOttll2010l) with pipeline version 2.4. 

The observations presented here are full spectral scans of 
Bands 6b and 7b, meaning they cover a frequency range of 
1573.4 - 1702.8 GHz (176.2 - 190.7 yum) and 1788.4 - 1906.8 
GHz (157.3 - 167.7 //m), respectively. These spectral scans con- 
sist of double-sideband (DSB) spectra with a redundancy of 4, 
which are deconvolved in to single-sideband (SSB) spectra . This 
procedure is outlined in iBergin et al.l ( 12010. this volumg) . We 
applied the standard HIFI deconvolution using the doDeconvo- 
lution task within HIPE with no channel weighting or gain cor- 
rection. Strong spurs and noisy DSB data sets were not included 
in the deconvolution and no fringing correction was applied. All 
data presented in this Letter are deconvolved SSB spectra. 

After the deconvolution was performed, the data were ex- 
ported to FITS format and all subsequent data reduction and 
analysis was performed using the IRAM GILDAS package. 
Main beam efficiencies for bands 6b and 7b were assumed to 
be 0.64 and 0.63, respectively. We estimate the typical RMS in 
both bands to be T^ x 0.9 K at the original spectral resolution. 



3. Results 

The SSB spectra for Bands 6b and 7b are given in Fig. [T] 
smoothed to a velocity resolution of ~4.5 km s ' and cor- 
rected for a WisR - 9 km s"' with the most prominent lines 
(peak Tmb S 7 K) labelled. Polynomial baselines of order 2 
are also subtracted from each spectrum. We find that these ob- 
servations are dominated by strong lines of CO, H2O, and OH 



as was reported by iLerate et aP (l2006h . With the higher spec- 
tral resolution of HIFI, we also detect additional strong lines 
of CH3OH, H2S, HCN, and HDO. Line identifications were 
made with the a id of the XCLASS prog rarrQ which accesses 
both the CDMS ('Muller et al.l l200lLl2005L http://www.cdms.de) 
and JPL (Pickett et al. 1998, http://spec.jpl.nasa.gov) molecu- 
lar databases. We list these transitions along with their inte- 
grated intensities in Table [T] Line intensities were measured 
using the CLASS data reduction and analysis software pack- 
age. In instances where there were blends, Gaussian profiles 
were fit to the lines and the results from the fitted profiles are 
reported; otherwise the total intensity is measured directly us- 
ing the BASE command. All line intensities were measured us- 
ing spectra smoothed to a velocity resolution of ~1 km s"'. 
Uncertainties in the integrated intensities, ct;, were computed us- 
ing the relation cr/(K km s"') = '\fN{6v)RMS where 6v is the res- 
olution in velocity space, N is the number of channels over which 
the intensity is measured, and RMS is the root mean square de- 
viation in the vicinity of the line. In addition to the lines listed 
in Table [1] we also detect many additional weak transitions of 
CH3OH, SO2, H2S, and H2O along with their isotopologues. 
Examples of several weak lines detected in Bands 6b and 7b are 
plotted in Fig.|2] Integrated line intensities for these weaker tran- 
sitions along with peak intensities for all lines will be reported 
in a later study. 

When comparing these spectra to other lower frequency 
HIFI bands, it is readily apparent that the line density is sig- 
nificantly diminis hed when compared to the low e r frequency 
bands (see e .g. iBergin et al] 12010. this volumet I Wang et al] 
I2OIO. inpressh . We estimate that the total fraction of channels 
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Table 1. Strong Lines in Band 6b and 7b 



1.85 10 1.854 10 1.8545 10 1.B55 10 1.B555 10 

Frequency (MHz) 



Fig. 2. A small sample of weak lines in Bands 6b (top panels) 
and 7b (bottom panels). All spectra are smoothed to a resolution 
of ~1 km s '. 



taken up by lines is ~0.23 in the lower frequency bands com- 
pared to ~0.07 in Bands 6 and 7. We reached these estimates by 
counting the number of channels in emission in the frequency 
ranges 858.1 - 958.1 GHz (Band 3b) and 1788.4-1898.5 GHz 
(Band 7b). We adopt these line density estimates as being rep- 
resentative of the low and high frequency bands, respectively. 
Although not formally presented in this Letter, a full spectral 
scan of Orion KL taken in Band 3b was also obtained as part of 
the HEXOS key program and used to estimate the line density 
here. These data were reduced in the same way as Bands 6b/7b. 
Both bands were smoothed to a velocity resolution of ~1 km s ' 
and any channel that had a value Tmb > 2.5 K (after baseline 
subtraction) was flagged as being in emission in 7b. This thresh- 
old is approximately what we have estimated as 3x the RMS 
in Tmb in Band 7b at a resolution of 1 km s"' (RMS ~ 0.8 K). 
Because the beam size, 6, decreases as a function of frequency 
(03fe~24"), fliis value was scaled to an equivalent RMS in Band 
3b using the following relation. 



' RMS^ 



TrmsJ^J 



(1) 



which assumes that the source is significantly smaller than both 
beam sizes. Thus the reduced beam size should be more coupled 
to the smaller spatial components (e.g. the hot core). One might 
therefore naively expect the line density to increase at THz fre- 
quencies. The opposite trend, however, is observed. 

4. Discussion 

One of the primary reasons for the reduced line density in the 
high frequency bands is the fall off" in emission from complex or- 
ganic molecules - in particular the "weeds" such as CH3OCH3, 
SO2, C2H5CN, and, of course, CH3OH. In Fig. [3] we present 
the number of emissive lines for select "weeds" as a function 
of frequency. To estimate these numbers we assumed LTE and 
predicted the emission for each species assuming T = 150K. We 
use the total colum n estimated for each molecule on the basis of 
IComito et alJ (1200 5) and in addition assumed a velocity width 
of 5 km s"^ . If the predicted emission was above O.IK then we 
counted the line as potentially emissive in our 100 GHz bins. In 
this fashion we counted Noak, which is shown in the figure. As 



Molecule 



Frequency 
(MHz) 



Transition 



j TMsdv Notes 
(Kkms-') 



Band 6b 



H20 


1574232.073 


64,3 -7 1.6 


70.8 + 6.7 




CH30H 


1586012.991 


85,1 - 74,1 


85.3 ± 3.0 


3 




1586013.008 


85,0 - 74,0 




3 


H2S 


1592669.425 


72,5 -7 1,6 


69.5 ± 3.4 




HCN 


1593341.504 


18-17 


1384 + 5.5 


1 


Ui^^S 


1595984323 


42,3 - 3l,2 


76.3 ±5.1 




CH3OH 


1597947.024 


96,0 - 85,0 


844 + 4.1 


3 




1597947.024 


96,1 - 85,1 




3 


H,S 


1599752.748 


42,3 -3l,2 


258.0 ± 6.4 




H,0 


1602219.182 


4l,3-4o,4 


959.0 ± 8.5 




H.34S 


1605957.883 


61,5 - 60,6 


116.0 + 6.0 


3 


H2'*'0 


1605962.460 


4l,3-4o.4 




3 


HoS 


1608602.794 


625 - 6| 6 


91.8 + 4.2 




CO 


1611793.518 


14-13 


4653.0 + 12.0 


2 


H2'*'0 


1633483.600 


22,1 -2i.2 


192.0 + 9.3 


3 


CH3OH 


1633493.496 


134,0 - 123,0 




3 


H2S 


1648712.816 


42,2-33,1 


124.9 + 4.9 


1 


"CO 


1650767.302 


15-14 


287.0 + 8.0 


3 


CH3OH 


1650817.827 


214,0-203,0 




3 


H.'^'O 


1655867.627 


2l,2 - lo.i 


75.6 + 8.2 


2 


H,0 


1661007.637 


22,1 -2l,2 


1008.0 + 9.2 


2 


H2"0 


1662464.387 


2l,2 - lo,i 


155.0 + 7.4 


2 


HiO 


1669904.775 


2l,2 - lo,i 


2266.0 + 10.4 


2 


HCN 


1681615.473 


19-18 


136.0 + 5.8 




CH3OH 


1682556.723 


105,1-94,1 


76.6 + 3.9 


3 




1682556.856 


105,0-94,0 




3 


HDO 


1684605.824 


61,5 - 60,6 


71.9 + 4.2 




Band 7b 


H2O 


1794788.953 


62,4 - 61,5 


969.0 + 8.5 




H.O 


1797158.762 


73,4 - 72,5 


648.0 + 49 




NH3 


1808935.550 


3i,i -2i,o 


158.0 + 6.4 


2 


NH3 


1810377.792 


32,1-22,0 


50.1 +8.3 


2 


H2'«0 


1815853411 


53,2-52.3 


91.0 + 5.5 




CH3OH 


1817752.285 


77,0 - 65,0 


35.7 + 3.2 




OH 


1834747.350 


2n,/23/2- - 1/2+ 


627.0 + 9.7 


2,4 


OH 


1837816.820 


2ni/23/2+ - 1/2- 


640.0+ 11.5 


2,4 


CO 


1841345.506 


16-15 


3820.0 + 10.4 


2 


H,S 


1846768.559 


61,6 - 5o,5 


237.0 + 5.8 




H2S 


1852685.693 


5l,4-42,3 


181.0 + 48 




H2S 


1862435.697 


52,4 -4,.3 


847 + 5.2 




H2S 


1865620.670 


33,0 - 22,1 


173.6 + 7.8 




H2O 


1867748.594 


53,2-52,3 


864.0 + 8.8 




H2O 


1880752.750 


63,4 - 7g,7 


135.0 + 54 




H2S 


1882773.396 


83,6 - 82,7 


37.9 + 4.0 




H,0 


1893686.801 


33,1 -4o,4 


265.0 + 5.9 




H2'**0 


1894323.823 


32,2 -3l,3 


161.0 + 6.4 




H2S 


1900140.572 


7 1,6 -7o,7 


83.80 + 6.7 


3 




1900177.906 


72,6 -7 1.7 




3 



Table 2. Notes: 1 - Line intensities were fit using Gaussian pro- 
files because of a blend. 2 - Lines contained self absorption. 3 - 
Subsequent lines with this note were severely blended and could 
not be separated by line fitting. The same integrated intensity is 
reported for both entries. 4 - Other OH transitions contributed 
to the integrated intensity reported. The strongest transition is 
given. 



can be seen, there is a general decrease in emission for all species 
but its particularly evident for CH3OCH3 and C2H5CN. CH3OH 
has a small factor of 2 decrease in the number of lines and, at the 
zeroth level, this is seen in our data which has numerous weak 
methanol transitions scattered throughout the band. 
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Fig. 3. Predicted number of lines with peak emission > 0. 1 K 
based on an LTE model in 100 GHz bins for select "weeds" as a 
function of frequency. 

Another possibility is that the dust emission from the hot 
core is optically thick in the high frequency bands; thus the dust 
would absorb all of the photons emitted from the molecules in 
the hot core. To explore this more closely we can examine the 
dust opacity expected within the hot core itself. Plume et al. 
(2010, in preparation) used multiple transitions of C"^0 and 
spectrally isolated the hot core. They estimate an N(C'**0) 



,23, 



1 .7x 10 cm '■ which yields a total H2 col umn of 3.4x 10 cm 



,-2 



assuming n(C"*0)/n(H2) ^ 1.7xlO-MFrerking et al.1 [T98l . 
Using the relation given in iHildebrandl (1 1 983l Equation 10), we 
estimate a t~0.1 at 171 ycca, putting it slightly lower than being 
optically thick. 

It is clear that there are other emission components in this re- 
gion as we see widespread emission from a variety of molecules 
in the high frequency bands. However, we still observe many 
molecules (CH3OH, H2O, HDO, and HCN) that have velocity 
components in their spectral profiles that are coincident to those 
expected from the hot core and other components (e.g. the out- 
flows). If these emission components do arise in the the hot core, 
it is likely that the molecular emission region must lie in front of 
any optically thick core. Given the presence of st rong physical 
gradients in the density and temperature profiles dWright et alj 
119961 lBlakeetal]ll996h and the fact that the dust is marginally 
optically thick, this is not unrealistic. 

A final contributor to the decrease in the line emission could 
be non-LTE excitation. At high frequencies there are a larger 
number of high excitation lines which could be more difficult to 
excite even at densities of 10^ - 10^ cm"-*. This needs to be more 
directly calculated using a molecule such as CH3OH with colli- 
sion rates that extend to temperatures greater than 200 - 300 K. 

5. Summary 

We have characterized the high frequency spectrum of Orion 
KL. We find that the spectrum is dominated by strong lines of 
CO, H2O, HDO, OH, CH3OH, H2S, HCN, and NH3. We also de- 
tect many weaker transitions of CH3OH, H2O, HDO, and SO2. 
We find that the line density is diminished in the high frequency 
bands when compared to the lower frequency bands and provide 
a number of explanations as to why this may be. 
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